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C/3 . Abstract. A new method is proposed to reduce the frequency noise of a local 

^~>' oscillator to the level of white phase noise by maintaining (not destroying by 

projective measurement) the coherence of the ensemble pseudo-spin of atoms over 

many measurement cycles. This method, which we call "atomic phase lock (APL)," 

^^ uses weak measurement to monitor the phase in the Ramsey method and repeat the 

^ . cycle without initialization of the phase. APL will achieve white phase noise as long 

^sO ! as the noise accumulated during dead time and the decoherence are smaller than 

1^ ' the measurement noise. A numerical simulation confirmed that with APL, Allan 

deviation is averaged down at a maximum rate that is proportional to the inverse of the 
^^ total measurement time, t~^. In contrast, current atomic clocks that use projection 

f— ^ ■ measurement suppress the noise only to the white frequency noise level, in which case 

the Allan deviation scales as r~^/^. Faraday rotation is one way to achieve weak 
measurement for APL. The strength of Faraday rotation with ^^^Yb^ ions trapped in 
a linear rf-trap is evaluated, and the performance of APL is discussed. The main source 
of the decoherence is a spontaneous emission, induced by the probe beam for Faraday 
''^ • rotation measurement. The Faraday rotation measurement can be repeated until the 

5t 1 decoherence becomes comparable to the signal-to- noise ratio of the measurement. The 

number of cycles for a realistic experimental parameter is estimated to be ^100. 



1. Introduction 

Many applications and experiments use the electromagnetic (EM) field to manipulate 
the quantum state of two- level systems (TLS). Improving the frequency stability of the 
EM field and TLS are of great importance because the ability to precisely and coherently 
control the population and phase of TLS plays a crucial role in many applications, 
such as nuclear magnetic resonance (NMR) spectroscopy and imaging, atomic clocks, 
magnetometers, quantum computers, and quantum simulators. The frequency (or 
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phase) of the EM field is usually referenced to the oscillator called the local oscillator 
(LO). 

A no-feedback approach, known as spin echo, has been developed in NMR to 
suppress the dephasing error [1]. The spin echo uses the phase of the LO as a reference 
and averages out the phase error of TLS by inserting vr-pulses. A similar technique 
can also be applied to isolate non-classical effects, such as entanglement with the noise 
environment, and is called "djTiamical decoupling" [2]. No-feedback approaches are very 
useful and easy to implement because the phase difference between LO and TLS does 
not need to be monitored. However, it does not suit applications such as atomic clocks 
and magnetometers, whose goal is to stabilize LO by using TLS as a reference. Spin 
echo also cannot maintain the phase error in the long term. Since we are pursuing the 
long-term stability of LO, we chose to use the feedback approach. Different methods 
are used depending on the kind of target and reference oscillators and whether the noise 
is suppressed to white frequency noise or to white phase noise. A summary of various 
feedback methods is shown in Table [H 



Target 


Reference 


White frequency noise 

ay oc r-V2 


White phase noise 

ay (X T^^ 


Laser 


Laser(LO) 


Transfer cavity 


Laser phase lock 


TLS 


LO 


NMR lock [3] 


Coherent magnetometry [1] 


LO 


TLS 


(Conventional) atomic clock 


"Atomic phase lock" 



Table 1. Various ways to match the frequencies of oscihators. t is the total 
measurement time, and ay is the Allan deviation (deviation of the target frequency 
from the reference frequency) of the normalized frequency (explained in Section [3]) . 
When both the target and reference oscillators are lasers, white frequency noise is 
often achieved by locking the target laser to the reference laser by matching the 
frequency of the resonant peaks. For example, the transfer cavity length is locked 
to the reference laser frequency, and then the target laser frequency is locked to one 
of the transfer cavity's resonances. The NMR lock looks at the NMR signal of the 
deuterated solvents, and the magnetic field is feedback controlled to keep the atomic 
spin resonant frequency constant [3] . An atomic clock is a typical example of an LO 
frequency being locked to the atomic spin resonance frequency. In both cases, noise 
was suppressed to white frequency noise, and (Jy was reduced at a rate of r~^" as 
derived in Section 12.11 



Conventional atomic clocks use the Ramsey method with projection measurement. 
Although this Ramsey method compares the phases of LO and TLS, it cannot achieve 
white phase noise over many measurement cycles because the phase of the atomic spin 
is destroyed and initialized at each cycle due to projection measurement. As these 
measurement cycles are repeated, the measurement noise accumulates at each cycle and 
the phase uncertainty grows as S^ oc r^/^. As a result, frequency stability decreases at 
a rate of ay oc r~^/^, which is characteristic of white frequency noise. To achieve white 
phase noise, we need to monitor the phase of the atomic spin over many cycles without 
destroying it. 
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We propose an experimental method called "atomic phase lock (APL)" to achieve 
white phase noise. This method combines weak measurements O |6] with the Ramsey 
method to monitor the phase difference while having the least effect on the coherence 
of the spin over many cycles. 

With APL, ay can be reduced at a faster rate up to r^^ as noise is suppressed to 
white phase noise. This ay oc r~^ is achieved when the phase of the target oscillator is 
locked to the phase of the reference oscillator. Although white phase noise is routinely 
observed when locking a target laser phase to a reference laser phase, feedback control 
of phase matching between TLS and LO has not been achieved. This is because TLS 
is a passive oscillator and monitoring the phase of TLS is difficult without affecting 
(destroying) the phase itself. 

This paper is organized as follows. In Section 2, atomic clocks that run with 
projection measurement are reviewed, and APL is introduced. In Section 3, we estimate 
the Allan variance for atomic clocks with projection measurement and with APL. 
In Section 4, we propose the ^^^Yb"^ ion trap with microwave transition as a clock 
transition for demonstration of the proof-of-principle. Then, the signal strength of 
Faraday rotation is estimated. The decoherence rate is calculated to estimate how many 
cycles of Faraday rotation measurement can be performed. In Section 5, we compare 
the stability of different types of atomic clocks using numerical simulation. In Section 
6, we discuss systematic shifts of the ^'''^Yb''' microwave clock, a comparison between 
APL and the spin-squeezed Ramsey clock, and the extension of APL to the optical clock 
transition. 

2. Ramsey method and APL 

2.1. Review of Ramsey method with projection noise 

The phase of the spin (relative to the LO phase) cannot be measured directly. The 
Ramsey double pulse method, in essence, measures the spin phase by mapping the phase 
information to the population ratio. Projection measurement of the population ratio 
gives us the phase. Note that the Ramsey method is often interpreted as a measurement 
of the average frequency shift during the time between double pulses. We instead 
view the Ramsey method as a measurement of the phase accumulated during the time 
between two pulses. This is because viewing it as a phase measurement is important for 
understanding the advantage of our proposed APL. If the Ramsey signal is expressed 
as a function of frequency, (p = J Auit) dt can be used to change the dependence on cj). 

Conventionally, the phase information represented by the population ratio has 
been measured via projection measurement. The Ramsey method with projection 
measurement is called "projection Ramsey" in this paper. This is to clarify the difference 
to APL, whose signal detection is via weak measurement. 

In this section, we will review the projection Ramsey method using Bloch 
representation [7]. The dynamics of a two-level system interacting with an EM field 
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are the same as the dynamics of a spin-| particle in a magnetic field [S], and we will 
only use this "spin" picture in the rest of this paper. The Bloch representation is 
useful for discussing the time evolution of the spin interacting with the field. For a 
more detailed explanation of the Bloch sphere in the context of atomic clocks, see, for 
example, Ref.[9]. 

The projection Ramsey sequence is shown in Figure [H The Cartesian coordinates 
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Figure 1. Step-by-step description of the Ramsey method with Bloch sphere pictures. 
The red arrow is the spin vector, and the blue arrow is the torque vector that rotates 
the spin. The angle in degrees is the phase of the EM field, and the angle in radians 
is the rotation angle of the spin by the EM field. Tc is the cycle time. 



of the Bloch sphere are designated as the u, v, and w axes. We limit our argument to 
the case where the phase of the spin is coherent (pure state). The ground and excited 
state correspond to the south and north poles respectively, and the phase of the spin 
corresponds to longitude. We use the frame that rotates at ujlo and the rotating-wave 
approximation. Then we consider the following two types of rotation. A spin is rotated 
around a torque vector that lies on the u — v plane when the resonant EM field is applied. 
A spin is also rotated along the longitudinal direction at a rate given by Aw, which is 
the frequency difference of LO and atoms. For clarity, we will express the phase of the 
EM field in degrees and the rotation angle in radians. 

Ramsey measurement proceeds as follows. (1) Repump all the atoms to \g) to 
correspond to aligning the spin to the coherent spin state pointing along —w. (2) Apply 
the strong resonant EM field to rotate the spin around -v by 7r/2 to point u. We define 
this phase of the EM field (LO) as 0° and assume the power is strong enough that the 
time taken for this rotation is negligible. (3) Wait for free precession time Tpp, and the 
frequency difference between LO and the atom rotates the spin around w. The angle 
between u and the spin corresponds to the phase difference of LO and atomic spin that 
is accumulated during Tpp. (4) Apply the second EM field with 90° phase shift to rotate 
the spin around u by 7r/2. Now the angle (f) is represented by the value along w;, which 
is a population ratio of the superposition state. Again, we assume that a negligible 
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Qm 




Figure 2. Graph of the mean detected signal as a function of the phase difference 
accumulated during Tpp in the Ramsey method. The solid and dotted lines show cases 
of +90° and —90° phase shifts respectively in step (4) in Figure [1] 



amount of time is taken for this step. (5) Projection measurement signal Q gives the 
measure of the population ratio between \g) and |e). The mean detected signal of a Q 
is shown in Figure [2j Signal Q{4>) is a sinusoidal function of the phase (j) and is given by 

Qr 



Q(0) 



:i + sin(0)) 



in Figure [2] as an error signal to 



The passive atomic clock uses this slope near q 
lock the frequency. 

The projection measurement in step (5) is normally done with electron shelving, 
first proposed by Dehmelt [ini [H] . Applying a laser beam appropriately polarized and 
tuned to a transition will scatter many photons if the atom is in one of the two level's 
states (a "cycling" transition) but will scatter no photons if the atom is in the other 
state. The population ratio can be measured by counting the number of bright ions. The 
value of the population ratio has a fundamental fluctuation of the quantum projection 
noise (QPN) [12]. QPN originates from the randomness in quantum projection and 
follows the binomial distribution (for the mathematical expression, see Eq. (]B.5|) ). This 
QPN is then reflected in the measurement of (p through Eq. ([T]). 

We shift the phase of the EM field by 90° to obtain the greatest sensitivity (largest 
slope) at around = 0. Applying two EM field pulses with the same phase is common, 
and this minor modification of the phase itself will not affect the performance of the 
atomic clock. 

The atomic clock's performance is normally improved by extending the free 
precession time (Tpp) of the Ramsey method, but Tpp is limited by the stability of 
the LO. In other words, the Tpp can be longer only as long as the phase difference 
accumulated in Tpp is guaranteed to be within ±7r/2. In the aim for a longer Tpp 
time, the stability of the laser (LO) has been improved by locking the laser frequency 
to a high-finesse cavity with ultra-low-expansion (ULE) glass as a spacer. However, the 
thermal noise of the mirror coating [13] is a hard limit to break through with the present 
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technology, and improving the Tpp by an order of magnitude would be difficult. In such 
a situation, APL could provide the alternate path to break through this limit to achieve 
a longer Tpp because APL is equivalent to lengthening the Tpp, as long as the atom 
coherence is maintained. 

2.2. Atomic phase lock 

The APL method we propose is structured as follows (Figure [3]). (1) Initialize the spins 
along —w. (2) Rotate the spin by a 7i/2 pulse with 0° phase. (3) Wait for the free 
precession time Tpp. (4) Rotate the spin by a 7i/2 pulse with 90° phase. (5) Estimate 
the +w component via the weak measurement. (6) Rotate back the spin by a 7r/2 pulse 
with -90° phase. (7) Repeat steps (3) to (6). 
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Figure 3. Stcp-by-stcp description of APL with Blocli splierc pictures. 



The differences from the projection Ramsey method are that the projection 
measurement is replaced by weak measurement and that, after the weak measurement, 
the spin is rotated back to where it was before step (4) and the process returns to step 
(3) without initialization. The point of APL is to monitor the phase without destroying 
the coherence of the spin so that the measurement noise does not accumulate after 
many cycles. A dispersion measurement serves as a weak measurement, and there are 
two ways to perform it: Mach-Zender interferometry [13] and Faraday rotation |15]. In 
Section HI we will discuss the Faraday rotation on spins of trapped ions. 



3. Allan deviation and stability of atomic clocks 

In this section, we will estimate the frequency stability of atomic clocks. For this, we 
use the Allan deviation, a widely used measure for evaluating the frequency noise of 
the LO, which is in general non- stationary and correlated. We first define terms and 
symbols. Fractional frequency deviation is defined as 

Au{t) v{t) - uq 



y{t) 



1^0 



^0 



(2) 
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where z/q is the clock frequency of atoms and assumed to be a constant. The Allan 
deviation of the fractional frequency cr,y(r) is defined as 



1 " _ 

Here, y means the time average of y over time r. The Allan deviation is the deviation 
from the previous data rather than from the mean. 

Now, we will estimate the Allan deviation of the atomic clock with the Ramsey 
method, using Figure [2l We follow the basic argument in Ref. [9] with a modification 
to view the Ramsey method as phase measurement. When there is a measurement 
uncertainty a^ , the corresponding uncertainty in </>, a'^ is linked by the slope of the 
spectrum as follows: 

d(t)~ a'*'' 
From Eq. ([1]), the slope at </> = is 



(4) 



d(p 2 

and Eq. (jl]) is given as 



(5) 



where SNR= Qmaxl^y^ is the signal-to- noise ratio of a single measurement. 

We define 0(r) as the total phase difference between LO and spin after time r. 
In contrast, A0 represents the phase difference accumulated during the free precession 
time of the projection Ramsey method. Allan deviation o"y(r) (r is total measurement 
time) is linked to the RMS deviation ^^(r) = y ([0('r) — (^('r))]^) as follows: 

The r dependence of 5^(r) is characterized by the type of noise (i.e., white frequency 
noise, white phase noise, etc.). For example, the repetition of the projection Ramsey 
measurement and feedback results in random walk phase noise because the spin phase 
is randomized by decoherence due to the projection measurement. Random walk phase 
noise is also called white frequency noise, and its deviation grows as ^^(r) oc a/t. This 
growth of the noise is due to "resetting (projection and initialization)" of the spin phase 
at each Ramsey cycle. In this case, 0(r) can be estimated as a sum of the phase 
measurement at cycles, 

</>(r) = ^ A</.„ (8) 

where iVc is the number of cycles given by Nc ~ t/^c and T^ is the time taken for a 
single Ramsey measurement cycle. Since A0j are non-correlated due to resetting of the 
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phase at each cycle, the ensemble average of the deviation (5^(r) is expressed using Eq. 
(ED as 



^y(^) = _,._o.TT. W7F- (10) 




S,{r) = v/([0(r)-(0(r))P) = ^ Mt)]^) = a^y^ ~ g^^^- (9) 

Here we consider only the case where r > Tc, and we used (0(r)) = because the LO 
frequency does not drift when it is stabilized to an atom. This growth of the deviation 
is analogous to the case of measuring a 400 m race track with a 10 cm scale stick. Every 
time 10 cm is measured and added to the total length, noise is added due to the finite 
width of the marker line. Finally, we get the Allan deviation of the atomic clock with 
the Ramsey cycle as 

1 
TTi/orSNR 

= ^— — d-. (Ill 

where Au = l/Tpp ~ l/Tc. Equation (fTO!) is the same as in Ref. [9], and we call this 
the "white frequency noise limit." When the SNR is limited by the quantum projection 
noise SNRqptv = 1/v^A^ {Na is the number of atoms), Eq. ( jTOj) is sometimes called the 
"QPN limit." 

From Eq. ( ITOj) . we can see that cry(r) decreases as I/a/T^;. Thus, it is better to use 
a longer free precession time, but normally, the length of Tpp is limited by the noise of 
LO. In other words, when Tpp is too long, the phase difference after Tpp becomes larger 
than ±|, and whether the signal is, for example, l/Air or 3/47r cannot be distinguished 
(see Figure E]). Therefore, in the regular Ramsey sequence, Tpp is kept sufficiently small 
that the phase signal is guaranteed to be within ±|. 

APL aims to overcome this limitation due to the noise of LO by use of the 
weak measurement. In essence, APL "peeks" at the phase with the least amount of 
decoherence, and feedback will keep the phase within ±|. As long as the coherence of 
the spin is sufficiently maintained, the APL achieves a single Ramsey measurement with 
long Tc, that is equal to the total measurement time r. As a result, Tc = r is achieved 
and Eq. ( ITOl) now becomes. 

We call this ideal case of the APL as the "white phase noise limit" in this paper. 

4. Possible experimental setup to achieve atomic phase lock 

4-1. Experimental setup 

We will first specify an experimental setup for a quantitative discussion. We chose a 
microwave atomic clock, using ^^^Yb"*" ions trapped in a linear rf-trap, as a possible 
proof-of-principle experiment. We will see the white phase noise limit until decoherence 
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degrades the performance of APL. We chose an ion trap for its long trapping time and 
chose the magnetic transition of the ground-state hyperfine sphtting as a clock transition 
for its long life time and long coherence time. For weak measurement, we will focus on 
Faraday rotation because it is somewhat easier to set up in an experiment. For simpler 
design and implementation of the experiment, we chose Ytterbium 171 ions that had 
the smallest nuclear spin of 1/2. 

An energy diagram of the Faraday rotation experiment setup is shown in 
Figure Hl^a). We assigned '^Si/2{F = 0,m/ = 0) as the spin-down state \l) and 
'^Si/2{F = 0,mf = 0) as the spin-up state |'|'). We use the transition between these two 
states as the clock transition. 




Oscillator 



b) Polarizer Linear rf-tfap 

(^^ — )) 

Probe 
beam 




Figure 4. (a) Energy diagram of ^^^Yb^, for the Faraday rotatfon measurenient. 
A probe beam for the coohng transition is shown. |e±) denotes the excited-state 
sublevels of dipole-ahowed optical transition, (b) Experimental setup for Faraday 
rotation measurement and the schematics of the feedback. 



We assume that one of the clock levels 1 1) has two dip ole- allowed optical transitions 
of 1 1) "H- |e±). This structure is common for the species of atomic clock including the 
optical frequency domain. Figure IDj^b) shows a setup for Faraday rotation measurement. 
The probe beam is linearly polarized and off-resonant with the detuning of A± for 
It) ^ |e±) transition, respectively. The probe beam uses the ^5*1/2 ^^ -P1/2 transition. 
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whose wavelength A is A = 369.5 nm, resonant atomic scattering cross-section ctq is 
(To = "2;^ = 6.5 X 10"^"^ m^, and natural line width F is F = 20 MHz [16]. In the following 
calculations, we assume that the cross-section of atom distribution in the probe beam 
Ap is Ap = 7r(200 /im)^ = 1.3 x 10"'^ m^ and the atom number Na is Na = 10^. For 
these parameters, on-resonance optical depth OD is 0D= Naao/Ap = 0.52. We set the 
Zeeman splitting 5 to be 5 = F. This detuning corresponds to applying 14.3 Gauss to 
trapped ions. 

We aim to measure the spin phase 0, which is related to w at step (5) of the APL, 
as 

(j) = arcsin w = arcsin(— — 1), (13) 

where A'^ is the number of atoms with spin up and is defined as N^ = ^{w+ 1). We 
need to be careful with the concept of N^. Each atom is in a superposition state, not 
in the eigenstate, so Nf means "the expected number of atoms with spin up if it were 
projected" and has the uncertainty given by quantum projection noise. 

4-2. Faraday rotation with atoms 

We will estimate the Faraday rotation angle and generally follow the notation in Ref. 
[T7] in this section. The phase shift for each circularly polarized mode {o'±) by a 1 1) 
atom can be written as 

^0 A±/F 

^^ = -J-pTTJ^Ajfy- ^'^^ 

Here we modeled the resonance with the Lorentzian spectrum. When ip^ and (f- are 
different, the polarization plane of a linearly polarized wave is rotated. This effect is 
known as Faraday rotation. The evolution operator of the Faraday rotation is expressed 



as a combination of the phase shift operators (see Appendix A ), 

^FR = exp(i(y9+n+) exp(2(y9_n_) 

= exp [i{ip+ + V5-)(n+ + n_)/2] t/pR, (15) 

where we define 

?7fr = exp(-27^(n+ -n_)/2). (16) 

and '-fz = ~{v+ ^ V-)- Since the common shift between a± modes does not affect the 
rotation angle, we neglect the factor between Up^ and t/pR and use Eq. flT6|l as the 
evolution operator of the Faraday rotation in the following. 
By using Eq. (TT4|) , 7^ becomes 

'yz = 2(3N^ = f3Na{l+w), (17) 

3 - ( A+/r A_/F \ ao 

^ " 1,1 + (2A+/F)2 1 + {2A^/Ty) 2Ap- ^ ^ 

Note that PNa is a value of 7^ when the atoms are equally populated to the | f) state 
and I i) state. For our projected experimental condition (see Section 14.41 ) of A = 0, 
i.e., A_|_ = — F and A_ = F, /3 is estimated to be /3 = —1.0 x 10"''' rad/atom. 
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For convenience in the following discussions, we introduce a general rotation 
operator of the polarization plane as 

R{e)=exp{i0{n+-n^)),, (19) 



where 6 is equal to the rotation angle of the linear polarization plane (see Appendix A ). 
?7fr = -R(— 7^/2), and the rotation angle of t/pR is — 7z/2. 

4-3. Signal of polarimeter 

The rotation angle is measured using a polarimeter that outputs the photon number 
difference between two orthogonal linearly polarized modes. When we set these modes 
h and v, the output of the polarimeter becomes 

Q = r]e{nh-n^), (20) 

where 77 is the quantum efficiency of the detector. Since the photons are randomly 
branched at the polarizing beam splitter, we need a finite photon number to distinguish 
in which direction and to what angle the polarization is rotated. To obtain the change 
from w = 0, the polarization angle of the incident pulse is set to it/ A + l3Na/2 from the 
h-axis. 

1^(0)) = R{7r/A + PNj2)\a)f,\0),, (21) 

where |a)/i is a coherent state of h-mode with a complex amplitude of a. Note that |ap 
is the mean number of photons in a probe beam pulse. 
The pulse after the transmission becomes 

mw)) = UF^mO)) = R (7r/4 - /3Naw/2) \a)h\0)v. (22) 

The mean value, the square mean value, and the variance of the error signal become 

{Q{w))= 7]e\a\^ siniP Naw), (23) 

{Q\w)) = r/2eV|2 (1 + |«|2 sm\f3Naw)) , (24) 

{AQ')= {Q\w))-{Qiw))' = rj'e'\a\'. (25) 

Since the average of the error signal is proportional to the w component for PNaW ^ 1, 
{Q{w)) can be used as an error signal. Eq. (l25l) is known as the shot noise. Note that 
we neglect the quantum projection noise in a collective spin system, in Eq. fl25l) . If the 
spin noise were concerned, w should be treated as a q-number and the variance would 
become (AQ^) = ?7^e^|Q;p(l + (3"^ N^\a\'^ {Aw"^)) . If the photon number is large enough 
to satisfy {Q{w))^ > (AQ^), we can recognize that the unbalance from Q = is caused 
by the atomic dephasing of w rather than the shot noise. For example, 

is required for the sensitivity of w; = 0.1 (see Figure [5l). We define the photon number 

A^p/i(SNR) as a number of photons that satisfies (Q(siIr))^ = (AQ^), so 

For our condition with A = 0, iV^ = 10^ and SNR= 10, A^ph(lO) = 9.3 x 10^ photons. 
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Figure 5. Shot noise y/{AQ^) and Faraday rotation signal {Q{w)). When {Q{w)) < 
^y{AQ'^), the signal is not distinguished from 0. 



4.4- Decoherence 

There are three types of decoherence mechanism: spontaneous emission, hght shift 
due to a probe beam, and backaction. The dominant decoherence is due to spontaneous 
emission and is discussed here. Decoherence due to hght shift and backaction is discussed 
in Appendix B 

For a hnearly polarized pulse with a complex amplitude a, {n±) = \aS^ j2. The 
total absorption probability by both circularly polarized modes for an atom due to a 
pulse illumination is 

CTqMp 



Pa 



%T 



i=± 



(2A±/r)2 



K) 



e\a\ 



where we set 



1 



1 



O-Q 



(2J 



(29) 



Vl + (2A+/r)2 l + (2A_/r)2y 2ylp■ 
For A = 0, £ = 1.0 X 10~^. 

Once an atom absorbs the photon, it will spontaneously emit the photon and is 
projected to 1 1)- This is the main source of decoherence with our experimental setup. In 
APL, we repeat the measurement with A'p/i(SNR) for iV^ep times. The total probability 
of decoherence is given from the total number of photons, as in Eq. fl28|) . and we obtain 

P,*°*'^'(SNR) = £Nph(SNR)N,,p. (30) 

Decoherence is safely permitted to be g^, so N^ep is given by 

1 



N, 



rep 



e A/p/,(SNR) SNR' 



Nj-ep as a function of A is shown in Figure O We can see that A', 



rep 



(31) 

is maximum at 

A = and has the value of 1.1 x 10^ when SNR= 10. This means that we can repeat 
APL for 100 cycles. 

Note that at A = 0, the probe beam is resonant to a ^6*1/2 (-^ 



l^mp = 0) -> 



Pi/2{F = l,rnp = 0) transition, but this transition is forbidden because AF = Anip 



0. 
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Figure 6. Nrep as a function of A. Na = 10^ SNR=10, T = 20 MHz, Ap = 1.3x 10"'' 
m^, and ctq = 6.5 x 10"^'*. 



5. Numerical simulation 

To confirm that the Allan deviation is reduced at rate r~^, we ran a simple numerical 
simulation. We first generated a free-run LO noise of the flicker frequency noise type. 
We used a noise level of cr^'"'^'^™"(r) = 1 x 10~^^ to model the noise of a crystal oscillator. 

The feedback schematics are shown in Figure H] (b). 

For the projection Ramsey method, we applied the following feedback scheme: 

A0« 



y projection 



(i-1) 
y projection 



Gf 



2'kvqTj 



FP 



(32) 



.(*) 



where yprojection is the updated normalized frequency at the end of the i-th cycle, Gp is 
the "proportional" gain, and A0^*) is the phase accumulated during the i-th cycle. We 
modeled the effect of projection measurement by resetting the phase at every cycle. 
For APL, feedback frequency is given by 



(i) _ (i-i) 
Vapl — Vapl 



G, 



(f)(.i) _ 0(i-l) 



(33) 



2Tfi'QTpp 

where 0*^*-* is the (total) phase measured via weak measurement. Note that this (/)'•*•' 
cannot be obtained unless the phase is maintained over (i) cycles. The measurement 
noise that corresponds to a'^ of the given SNR is included in A0 or 0. We assumed that 
the dead time was zero for both projection Ramsey and APL. For APL, we included 
the effect of the decoherence by resetting the phase every 100 cycles. 

The results of the numerical simulation are shown in Figure [3 This graph show 
that the Allan deviation is averaged down at a rate of oc r^^^^ for a feedback cycle with 
the projection Ramsey method. For APL, the stability decreases as r~^, from r = 0.5 
to 50 s. The phase is reset every 50 s (100 measurement cycles), and thus the stability 
decreases as r~^/^ from 50 s, which is the same slope as that of the projection Ramsey 
method. 

SNR=10 was chosen to minimize the decoherence rate during APL. For comparison, 
the SNR of the projection Ramsey method was also set to 10, representing a case when 
SNR is limited by technical noise [I2] at the 10 % level. 
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Figure 7. Numerical simulation of the Allan deviation for the projection Ramsey 
method, APL, and free-run LO. The phase of the spin is reset every 100 cycles (50 s). 
Step time width dt = 50 ms, carrier frequency vq = 12.6 GHz, SNR=10, cycle time 
Tc = Tpp = 500 ms, and proportional gain Gp = 1. 



6. Discussion 



For APL to achieve white phase noise, we need to trap a large number of ions in an 
optical path and also keep the dead time sufficiently small compared to Tpp. Dead time 
is the time taken by steps (2), (4), (5), and (6) in Figure El For the number of atoms, 
we chose Na = 10^ to have a decent OD. If we try to confine ions in Ap = 1.3 x 10~^ m^ 
and assume the inter-distance of ions to be 30 /xm, then the trap length needs to be at 
least 20 cm to have 10^ ions. We can avoid making a long trap by forming a bad cavity 
along the probe beam, around the ions of finesse ~100, and reduce the number of ions 
to 10"^ and the size of the trap to 2 mm. As for the dead time, its effect in a case with 
projection noise has been thoroughly investigated by Dick [18]. For APL, we need to at 
least keep the dead time low enough not to be limited by the so-called Dick limit. In 
other words, APL will work only as long as the noise accumulated during the dead time 
is negligible. A quantitative discussion would require careful evaluation and simulation, 
but this is beyond the scope of this paper. 

We have discussed the stability of APL based on a simulation assuming that atomic 
clock frequency does not change over time. When atomic clock frequency is unstable, 
the stability of the atomic clock will be limited by the stability of the atomic clock 
frequency itself. The systematic uncertainty of the ^^^ Yb"^ ions, cooled by helium buffer 
gas, is reported to be 6 mHz [19]. Even with this rather high systematic uncertainty, 
we can still expect the Allan deviation to decrease to cr^, = 5 x 10~^^ or lower without 
further cooling. 

The experimental setup for APL coincides with that of spin squeezing. Wineland 
proposed using the spin-squeezed state to overcome quantum projection noise [201 l2T] . 
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and spin squeezing for microwave clock transitions has been demonstrated [221 123] • The 
comparison between APL and spin squeezing is an interesting topic and will be published 
elsewhere. Here, we briefly mention that both use the weak (dispersion) measurement, 
but squeezing aims to improve the SNR in one measurement while APL initially uses a 
lower SNR but aims to improve the stability over many cycles by preserving the phase 
longer and avoid the QPN limit in the long term. 

In a similar spirit, to improve the (short term) stability of the atomic clock, an 
optical active clock has been proposed [211 ESj- However, the design of the optical 
active clock consists of an optical cavity, and long-term stability is again limited by the 
thermal noise of the mirror coating. 

Note that APL primarily contributes to improving the stability, and not directly to 
the accuracy of the atomic clock. However, improvement of the stability would enable 
the experimental estimation of the accuracy to be much shorter and indirectly contribute 
to improving the accuracy. 

Once APL is demonstrated with the microwave clock, the same principle should 
be applied to the neutral atom clock trapped in the optical lattice. We propose to use 
rf-ion traps for proof-of-principle demonstration with a microwave clock, but there will 
be difficulties in applying APL to such traps. First, buffer gas cooling will no longer 
be sufficient for ions to stay within the Lamb-Dicke limit because the wavelength is 
much shorter. Second, AC stark shift due to a trap rf-fleld will shift the optical clock 
frequency. An optical lattice clock is a candidate for avoiding those two problems. The 
time scale for applying APL will be shorter (up to ~1 sec) because the trapping time 
of the optical lattice clock is short (~lsec) compared to an rf-ion trap (~days). 
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Appendix A. Phase shift operator and rotation operator of polarization 
plane 

An evolution operator that gives the phase shift is exp{i(p±n±) , where n± is the number 
operator for a± mode. This can be seen by the following equation. After unitary 
conversions with the phase shift operators, the photon-annihilation operator for a± 
mode, say a±, becomes 

exp{—iip±n±)a±exp{iip±n±) = a±e*'''*, (A.l) 

where the commutation relation between the number and annihilation operator 
[n±, a±\ = —a±, and the Baker-Hausdorff lemma exp(Ai?)/lexp(— A-B) = A + X[B, A] + 
(A^/2!) [B, [B,A]] + ■ ■ ■ are used. Since the annihilation operator in quantum optics 
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corresponds to the amplitude of the EM field in classical electrodynamics, it can be said 
that the phase shift operator exp{iip±n±) gives the phase shift. 

An evolution operator that gives the rotation of polarization plane takes the form 
of Eq. (TT9|) . This can be seen by the following equation. We introduce two linear 
polarization modes h and v, which are related to the a± mode, as 

a^ = — ^. (A.2) 

After the unitary conversion with the rotation operator R{6), the photon-annihilation 
operator ah,v becomes 

R\e)ahR{e) = "+^' ^/"^ ' = ah cos 9 - a^ sin 9, (A.3) 

v2 

id —iO 

R\9)ayR{9) = '^ j3^ — = Oh sin 9 + a^ cos 6*, (A.4) 

V2i 



where Eq. (jA.ip is used. Thus, it can be said that the rotation angle of the polarization 
plane is 9. 

Appendix B. Light shift and backaction 

The probe beam affects the spin state in two ways: light shift (AC Zeeman shift) and 
imbalance fluctuation of the right and left circularly polarized component. We call the 
latter case "backaction." Both cases rotate the spin along the longitudinal u—v direction, 
but the difference is that light shift rotates the spin direction at a constant rate while 
backaction rotates randomly, resulting in increase of the deviation along the longitudinal 
direction. Light shift is zero for the case with A = because the contributions from 
'^Pi/2{F = 0, mp = ±1) cancel each other. In the following, we will show that the effect 
of backaction is also negligible compared to the decoherence due to absorption, which 
has already been discussed. 

To estimate the backaction, we need to treat the number of up spins (represented 
by w) as a q-number. For this purpose, we need to treat {u, v, w) as collective spin 
operators [26] . 

I Na 

« = ^ E (l5> t)^^{9, l\ + \9, lh{9, t I) (B.l) 

* =^Y.i-\9A)n{9,i\ + \9,i)^^{9A\) (B.2) 

I Na 

^ = TT E d^' t)u{9, t\-\9, l)u{9, i I) , (B.3) 

where \9,t)i (|fl')i)j) represents the state that the z-th atom is \9,t) (Ifi')i))- Their 
commutation relations are as follows: 

2i 2i 2i 

[u,v] = —w, [v,w] = —u, [w,u] = —v. (B.4) 
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From one of the commutation relations, we obtain an uncertainty relation 

{Av'){Aw') > ^. (B.5) 

a 

In the following discussion, we consider (u) = (n^) = 1. The minimum uncertainty state 
(Aw^) = (Aw^) = 1/Na corresponds to the quantum projection noise limit. 

After a probe pulse illumination, the atomic state operator evolves from A to A' as 

A' = {^{0)\Ul^AUFnm0)), (B.6) 

where 

Ufr = exp (-^/3iVa(l + w)in+ - n_)) (B.7) 

is the evolution operator of the Faraday rotation treating w as a q-number. By using 
the equations of motion under a fictitious magnetic field, 

UI^^uUfr = u cos (2/3(n+ — n_)) — v sin (2/3(n_|_ — n_)) , (B.8) 

UI^vUfk = u sin (2/3(n+ - n_)) + v cos (2/3(n+ - n^)) , (B.9) 

uIrU^Ufr = w, (B.IO) 

and the fluctuation of the circularly polarized component for linearly polarized coherent 
hght 

mO)\in+-n^)mO)) =0, (B.ll) 

(•<A(0)|K-n_)XO)) = |ap, (B.12) 

we obtain the increase of the variance 

(Am'2) = 4/3>|2({;2) (B.13) 

(A{)'2) = {Av^) (l - 4/32|ap) + Al3^\a\\ (B.14) 

In the second term of Eq. (]B.14p on the right side, 4/3^ |ap represents the backaction 
noise. For our experimental parameters, variance is increased along the longitudinal 
direction by 4/3^ |ap = 3.7 x 10~^*^ per measurement of SNR=10. The backaction is 
negligible under the condition of 4/3^|«piV(j <C 1. This condition is satisfied as long as 
the sensitivity of the Faraday rotation is well below the quantum projection noise limit. 
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